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SUMMARY 


The  objectives  of  this  study  were  to  review  the  use  of  vanes  for  sensing 
flow  direction  and  to  design  a  fleldable  vane  flow  direction  gage  for  use  in 
experiments  with  precursed  blast  waves,  and  to  recommend  a  test  program  for 
the  gage. 

Estimates  of  the  requirements  for  a  flow  direction  sensor  were  made  from 
the  results  of  an  existing  numerical  calculation  of  a  precursed  blast  wave. 
The  precursed  flow  field  Is  very  complex,  and  contains  two  large  vortices  and 
mechanisms  for  generating  smaller  scale  vortices  near  the  surface.  There  are 
high  speed  flows  (Mach  3)  present  In  a  Jet  near  the  ground,  and  the  response 
time  for  a  flow  direction  sensor  must  be  short,  0.2  ms  or  less  If  possible. 

A  rotating  wind  vane  sensor  was  examined  and  found  to  have  a  response 
time  that  was  too  long.  Alternative  means  of  sensing  flow  direction  were 
investigated  and  that  selected  for  a  gage  was  measuring  the  direction  of 
bending  of  a  cantilever  beam  protruding  into  the  flow  by  use  of  strain  gages. 

The  proposed  gage  has  four  cylindrical  hollow  cantilever  beams  in  a 
cruciform  configuration  on  the  end  of  a  supporting  shaft.  The  four  output 
signals  describe  the  flow  direction  in  the  vertical  plane  containing  the  gage 
and  ground  zero,  and  the  horizontal  plane.  Each  cylinder  has  a  natural 
frequency  of  about  8300  Hz.  Because  it  is  underdamped,  the  upper  usable 
frequency  is  about  2080  Hz,  with  a  rise  time  of  0.12  ms. 

The  test  program  recommended  includes  static  calibration,  tests  of 
response  versus  flow  angle  in  a  0.61  metre  diameter  shock  tube  at  a  shock 
overpressure  of  about  172  kPa  (25  psl),  supersonic  wind  tunnel  tests  at  Mach 
1.5,  2.0,  and  2.5,  and  tests  In  the  TRW  0.43  metre  shock  tube  at  higher  shock 
pressures  with  and  without  dust  in  the  flow. 
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SECTION  1 


INTRODUCTION 


The  blast  wave  from  a  nuclear  explosion  is  affected  by  the  hot  layer  of 
gas  produced  on  the  ground  surface  by  the  thermal  radiation  from  the  fireball. 
This  hot  layer  with  its  higher  sound  speed  alters  the  propagation  of  the  blast 
wave  so  that  a  wave  called  the  precursor  extends  into  the  hot  layer  in  advance 
of  the  main  shock  wave.  The  presence  of  the  hot  layer  causes  the  development 
of  a  very  complex  flow  field  with  large  and  fast  changes  in  flow  variables 
with  time  and  position.  For  targets  on  the  ground,  the  blast  loading  can  be 
radically  altered.  The  development  of  an  understanding  of  precursed  blast 
waves  and  a  capability  to  predict  the  flow  fields  is  essential  for  loading 
predictions  on  targets. 

Experimental  and  computational  studies  of  precursed  flow  show  that  large 
changes  occur  in  the  direction  of  flow.  At  this  time  there  is  no  suitable 
sensor  for  measuring  flow  direction  in  the  precursed  blast  wave  experiments. 
Experiments  have  been  conducted  and  more  are  planned  that  utilize  a  helium 
layer  at  the  ground  surface  to  simulate  the  high  sound  velocity  hot  layer. 
Flow  direction  sensors  are  needed  to  measure  flow  direction  versus  time  in  and 
above  the  helium  layer.  The  data  obtained  will  be  used  with  data  from  other 
types  of  sensors  for  diagnosis  of  the  blast  flow  field  and  for  comparison  with 
code  predictions. 

The  tasks  to  be  performed  in  this  study  were  as  follows: 

1.  Conduct  a  literature  survey  to  review  the  use  of  vanes  for  sensing 
flow  direction.  Determine  the  characteristics  of  the  precursor 
environment  to  understand  its  severity.  Establish  the  minimum 
requirements  for  a  vane  flow  direction  sensor. 

2.  Design  a  fieldable  flow  direction  sensor  which  is  ready  for 
construction  and  evaluation,  and  estimate  its  cost.  Recommend  a  program 
to  test  and  certify  the  sensor. 


1 


SECTION  2 


PRECURSOR  ENVIRONMENT 


Precursors  that  occurred  on  nuclear  explosions  produced  enhanced  loading 
on  ground  targets  and  Increased  the  radii  for  damage  to  some  types  of  military 
equipment  by  as  much  as  50  percent.  The  precursor  flow  fields  contained  dust, 
with  the  amount  of  dust  varying  with  test  location,  ground  surface,  and  yield. 
The  dust  added  to  the  loading  on  targets,  and  caused  damage  directly  by 
sandblasting  surfaces  and  plugging  gage  openings.  Dynamic  pressures  were 
increased  by  factors  ranging  from  two  to  four  above  those  that  would  occur  for 
undisturbed  blast  waves.  Dynamic  pressure  waveforms  were  characterized  by  a 
distorted  shape  and  large  and  rapid  oscillations.  The  environment  was  severe 
for  both  targets  and  gages. 

Recent  work  has  shown  that  the  primary  factor  affecting  precursor 
production  and  its  characteristics  is  the  high  sound  speed  layer  produced  by 
the  thermal  radiation  from  the  explosion.  Experiments  using  a  helium  layer  to 
produce  the  high  sound  speed  layer  on  a  surface  have  been  conducted  in  shock 
tubes  and  on  high  explosive  tests  where  charges  have  been  on  the  surface  and 
suspended  above  the  surface.  Precursors  were  produced  with  shock  wave 
configurations  and  pressure  waveforms  similar  to  those  observed  on  nuclear 
tests . 

Reference  1  describes  the  fluid  dynamics  of  precursors  and  the 
requirements  for  simulating  them  using  a  helium  layer.  The  results  of  an 
inviscid  2-D  numerical  calculation  of  a  clean  air  (no  dust)  blast  flow  field 
with  a  helium  layer  modeling  the  high  sound  speed  layer  are  presented  to 
display  the  characteristics  of  the  precursor  flow  field.  The  results  are 
presented  in  terms  of  a  1-KT  nuclear  explosion  at  a  height-of-burst  of  61  m 
(200  ft)  and  a  helium  layer  0.061  m  (0.2  ft)  thick. 

Figure  la  depicts  the  isopycnics.  Figure  lb  the  isobars.  Figure  2a  the 
velocities  fixed  at  the  moving  stagnation  point,  and  Figure  2b  the  velocities 
in  laboratory  coordinates  in  the  precursor  region  at  a  time  of  0.160  s/KT^/^. 
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These  figures  were  taken  from  Reference  1,  except  for  Figure  2a  (Ref.  2).  The 
incident  shock  has  run  forward  of  the  main  shock  in  the  high  sound  speed  layer 
(TL)  at  the  surface,  forming  the  oblique  precursor  shock  (P) .  The  precursor 
shock  P  causes  the  gases  in  the  high  sound  speed  layer  to  rise  upward,  which 
in  turn  induces  air  to  move  in  underneath  the  lofted  gas.  This  action  forms  a 
high  velocity  (Mach  3)  wall  jet,  which  is  very  obvious  in  Figures  2a  and  2b. 

The  explanation  of  the  figures  from  Reference  1  follows,  since  that 
reference  may  not  be  readily  available. 

"In  precursor-fixed  coordinates  (see  Figure  2A) ,  air  from  above  the 
thermal  layer  flows  through  the  oblique  precursor  shock  P  and  is  thus 
deflected  upward,  is  turned  approximately  parallel  to  the  ground  by 
rarefaction  wave  RW^,  and  is  turned  toward  the  ground  by  the  turning 
shock  T.  This  flow  then  partially  stagnates  at  the  'moving  stagnation 
point'  SP,  and  is  accelerated  forward  by  rarefaction  wave  RW2  which 
supports  the  wall  jet  flow.  In  effect,  rarefaction  waves  RW^  and  RW2 
correspond  to  the  pressure  fields  associated  with  two  large-scale  vortex 
structures,  Vj  and  V2,  with  positive  vorticity.  They  induce  counter¬ 
clockwise  rotational  velocities  which  add  to  the  background  motion. 
Hence,  velocities  in  laboratory  coordinates  are  increased  underneath  each 
vortex  and  decreased  above  each  vortex.  Velocities  on  the  wall  tend  to 
peak  and  pressures  drop  under  each  vortex.  The  above  features  correspond 
to  the  laminar  evolution  of  the  wall  jet. 

"Actually,  the  wall  jet  evolves  chaotically.  It  consists  of  a  free  shear 
layer  (FSL)  on  top  of  the  jet  and  a  wall  boundary  layer  which  was  not 
modeled  here.  The  free  shear  layer  is  hydro-dynamically  unstable  in  the 
Kelvin-Helmholtz  sense.  Small  perturbations  amplify  by  an  inviscid 
mechanism,*  and  agglomerate  in  coherent  vortex  structures  —  this  code's 
representation  of  the  turbulent  mixing  in  the  free  layer  region.  The 
wall  jet  acts  like  a  supersonic  channel  flow.  Area  variations  induced  by 
the  coherent  structures  on  FSL  create  'barrel  shocks'  and  rarefaction 
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waves  embedded  in  the  wall  Jet.  Corrugations  on  top  of  the  free  shear 
layer  induce  additional  oblique  shocks  in  the  flow  above  the  jet. 

"These  shocks  Induce  an  additional  nonlinear  wrenching  of  the  free  shear 
layer.  In  the  latter  stages,  the  mixing  in  the  free  shear  layer  extends 
all  the  way  to  the  ground.  This  spreads  the  jet  momentum  vertically,  and 
eventually  leads  to  the  turbulent  self-destruction  of  the  wall  jet  in  the 
precursor  cleanup  phase." 

The  discussion  above  indicates  that  there  are  at  least  two  large  scale 
vortices  Vj^  and  V2  to  be  expected  in  the  precursor  flow  field  and  smaller 
scale  vortices  generated  by  the  wall  jet  interaction  with  the  free  shear  layer 
above  it  and  the  boundary  layer  below  it.  Large  changes  in  flow  direction  may 
occur  for  a  sensor  placed  so  that  a  vortex  passes  over  it. 

The  flow  field  in  Figure  2b  displays  a  variety  of  flow  direction  changes. 
Along  the  precursor  shock  P  the  vertical  angles  of  flow  vary  from  about  40 
degrees  near  ground  level  to  about  60  degrees  near  the  main  shock.  The 
velocity  magnitudes  are  relatively  small  compared  to  the  wall  jet  velocities. 

At  182  m  (597  ft)  at  a  height  of  about  3.7  m  (12  ft)  in  the  vortex  Vj  the 
velocity  changes  direction  from  an  angle  of  about  15  degrees  with  respect  to 
the  horizontal  to  90  degrees  at  180  m  (590  ft)  to  180  degrees  at  178  m  (583 
ft),  a  change  of  165  degrees  in  4.27  m  (14  ft).  Assuming  a  velocity  relative 
to  the  ground  of  800  m/s  (Figure  2a),  the  change  in  angle  of  165  degrees  will 
occur  in  about  0.0053  s  as  the  flow  field  sweeps  past  a  sensor  located  at  180 
m  (590  ft),  corresponding  to  a  rate  of  3.1  x  10^  degrees  per  second. 

Lower  in  elevation  under  the  forward  edge  of  the  vortex  Vj  at  182  m  (597 
ft)  at  a  height  of  about  1.5  m  (5  ft)  the  velocity  angle  change  is  about  40 
degrees  over  a  distance  of  2.2  m  and  a  time  of  0.0027  s. 

Between  153  m  (503  ft)  and  155  m  (510  ft)  near  the  upper  part  of  the  wall 
jet  at  a  height  of  about  1.8  m  (6  ft)  the  velocity  changes  direction  about  30 
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degrees  over  a  distance  of  2.2  m.  Assuming  the  velocity  relative  to  the 
ground  Is  800  m/s,  the  time  for  the  change  Is  0.0027  s. 

In  the  wall  Jet  near  the  surface  the  angular  changes  In  direction  shown 
In  Figure  2b  are  about  10  degrees  or  less.  However,  the  numerical  model  does 
not  Include  a  boundary  layer  model  and  Is  limited  In  resolution,  so  that 
larger  angular  variations  may  appear  In  a  more  detailed  modeling. 

The  maximum  rate  of  change  noted  above  was  for  a  change  of  90  degrees  In 
0.00275  s  for  the  large  scale  vortex  Vj.  At  0.040  KT  or  20  tons  of  HE,  the 
time  Interval  would  be  0.00094  s,  and  for  0.001  KT  or  1000  lbs  of  HE,  0.000275 
s.  To  measure  a  change  of  90  degrees  In  0.00275  s  In  a  1-KT  yield  with  a 
resolution  of  5  degrees  would  require  a  sensor  response  time  of  0.15  ms. 

Considering  the  angular  rates  calculated  above  and  the  expectations  of 
additional  vortlclty  In  an  experiment,  the  response  time  for  a  vane  flow 
direction  sensor  should  be  about  0.2  ms  or  less,  "as  fast  as  possible." 

The  numbers  above  were  obtained  by  reading  from  the  figures.  Improved 
estimates  could  be  obtained  by  processing  the  hydrocode  data  to  read  out 
velocity  direction  and  magnitude  versus  elevation  at  several  ground  ranges. 
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SECTION  3 


ROTATION  WIND  VANE  RESPONSE  CHARACTERISTICS 

Simple  wind  vanes  have  been  In  use  for  many  years  to  Indicate  wind 
direction.  Such  vanes  consist  of  an  aerodynamic  surface  attached  to  an  arm 
that  Is  free  to  rotate  on  a  pivot,  and  a  mass  at  the  end  of  the  arm  away  from 
the  vane  that  statically  balances  the  vane.  A  properly  designed  vane  will 
align  itself  parallel  to  the  air  flow,  and  If  displaced,  will  return  to  its 
former  equilibrium  position  through  a  few  oscillations  of  rapidly  diminishing 
amplitudes.  The  design  of  such  wind  vanes  has  progressed.  Vanes  have  been 
used  in  a  variety  of  meteorological  and  aircraft  applications  for  a  wide  range 
of  flow  conditions. 

A  small  wind  vane  that  rotates  to  align  itself  parallel  to  the  air  flow 
was  examined  as  a  candidate  for  a  flow  direction  sensor.  The  response  of  such 
a  vane  for  small  deflections  will  be  used  to  judge  its  promise  for  blast 
measurements. 

For  small  deflections  the  equation  describing  the  vane  response  is 
analogous  to  that  for  a  spring-mass  system  with  viscous  damping  undergoing 
linear  free  oscillations.  The  derivation  of  the  equation  is  given  in  both 
References  3  and  4,  and  is  summarized  below. 

Figure  3  shows  the  schematic  of  a  simple  vane.  The  wind  vector  U  makes 
an  instantaneous  angle  6  with  respect  to  the  vane.  The  vane  experiences  an 
aerodynamic  force  F  acting  normal  to  its  surface  at  its  center  of  pressure,  a 
distance  r  from  the  axis  of  rotation.  The  vane  velocity  at  the  center  of 
pressure  is  rdS/dt.  In  coordinates  moving  with  the  vane,  the  air  velocity 
relative  to  the  vane  is  Uy,  and  the  instantaneous  angle  of  attack  By  is,  for 
small  S, 
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V 


o  +  I  M 
^  ^  U  dt 


(1) 


8 


When  the  vane  is  displaced  from  its  center  of  equilibrium,  the  normal 
force  F  produces  a  returning  torque  Fr.  For  small  B,  “  lUyl  ■  U.  Then 
the  aerodynamic  force  F  normal  to  the  vane  is 

F  =  (pU2/2)A  Cl  cos  By*  (2) 

where  Cl  Is  a  lift  coefficient,  A  is  the  projected  vane  area  in  the  plane  of 
zero  angle  of  attack,  and  p  is  the  density  of  the  air. 

For  small  Byf 


(3) 


Let  a 


Then  F  =  (pu2/2)A  a  By  (4) 


Under  dynamic  equilibrium. 


d^B 

dt^ 


=  -F  r 


(5) 


where  I  is  the  moment  of  inertia  of  the  oscillating  vane  about  the  rotation 
axis. 

A  unit  angle  torque  is  defined  as 


N  =  Fr/By  =  (pU2/2)  A  r  a 


(6) 


and  upon  its  substitution  in  equation  (5)  yields 


I  ^  .  .NBy  -  -N($  If) 

dt 


l^  +  ^^  +  NB*0 

dt^  “ 


(7) 
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The  aerodynamic  damping  constant  D  Is  given  by 


D  =  Nr/U 


Thus  the  equation  of  vane  motion  can  be  written  as  follows: 


d^B  +  D  dB 
dt2  I 


0 


(8) 


(9) 


In  terms  of  the  angular  natural  frequency  and  the  damping  ratio  the 
equation  Is 


2^0) 

n 


dt 


2 

+  0)  B 
n 


0 


(10) 


where 


0)^  =  (N/I)^  =  ([pU^/2]  Ara/I)^ 


(11) 


^  2 — r  =  (pr^Aa/8I)*^  (12) 
2(IN)^ 

The  natural  frequency  varies  linearly  with  wind  speed  U,  while  ^  is  a 
property  of  the  vane  alone,  providing  that  the  speed  range  is  small  enough 
that  the  lift  characteristics  do  not  change. 

The  time  lapse  for  oscillations  to  decrease  from  a  peak  amplitude  of  Bq 
at  t  =  0  to  some  specified  peak  amplitude  Bg  at  tj^  is  defined  as  the  vane 
response  time  (Ref.  4).  This  vane  response  time  can  be  obtained  from  the 
expression  for  the  amplitude  of  the  envelope  of  the  oscillations. 


and 


at  t  =  tR,  £n  (6^/e_)  =  Cw 

o  i\  n  K 

tn  (6^/8,) 

H 

Upon  substituting  the  relations  in  equations  (11)  and  (12)  for  and  c  > 
the  result  obtained  is 


(14) 


41  ln(6^/6j^) 
pAr^aU 


(21)^  In  (gjj/Bg) 
(pAra)^  CU 


(15) 


The  vane  response  time  increases  with  increasing  moment  of  inertia  I,  and 
decreases  with  air  velocity  U.  Fortunately  for  blast  wave  applications  the 
changes  in  flow  direction  are  expected  to  occur  in  the  early  portion  of  the 
blast  wave  where  the  air  velocities  are  near  maximum. 
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SECTION  4 


TIME  RESPONSE  FOR  A  CANDIDATE  VANE  DESIGN 


A  vane  for  use  In  precursed  blast  waves  will  be  subjected  to  supersonic, 
transonic,  and  subsonic  flows.  The  stability  of  a  vane  over  such  a  span  of 
flow  regimes  is  of  particular  concern.  The  candidate  vane  design  selected  was 
derived  from  Reference  5,  which  describes  tests  of  different  vane 
configurations  for  a  range  of  Mach  numbers  from  0.3  to  3.0  and  for  angles  of 
attack  ranging  from  -3  to  +24  degrees.  The  report  describes  a  vane  that  was 
dynamically  stable  over  the  entire  Mach  number  range.  Figure  4  shows  the 
details  of  the  vane  configuration  and  the  supporting  boom.  Two  sets  of  vanes 
are  shown,  one  to  measure  angle  of  attack  and  the  other,  sideslip.  The  vane 
configuration  shown  is  the  only  one  found  in  the  literature  search  that  was 
tested  over  such  a  range  of  Mach  numbers  and  angles  of  attack. 

For  application  to  blast  wave  measurements,  the  configuration  shown  in 
Figure  4  was  scaled  down  by  a  factor  of  about  seven  to  the  smallest  boom  that 
seemed  practical  for  use  in  the  field.  Figure  5  shows  the  configuration 
selected.  The  vanes  are  to  be  made  of  magnesium  alloy,  selected  for  low  mass 
and  high  stiffness,  the  shaft  of  hardened  steel,  and  the  counter-weight  of 
tungsten  to  minimize  size. 

Parameters  for  calculating  the  time  response  and  natural  frequency  are  as 
follows : 

I  =  5.75 

A  =  1.06 

r  =  1.1 

a  =  1.6 

The  parameter  a  was  calculated  using  the  relation 


X  10-8  kgm2 
X  10”^  m^ 

X  10-2  m 


a  =  2ir 


AR 

2  +  AR 


(16) 
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ALL  DIMENSIONS  IN  INCHES 

Figure  4.  Schematic  of  flight  path  accelerometer  with  movable 
vanes  from  Reference  5. 
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Figure  5.  Rotating  component  of  candidate  vane  sensor. 


where  Ar  is  the  aspect  ratio,  which  is  0.71  for  the  particular  vane 
configuration.  This  relation  was  found  to  be  in  fair  agreement  with 
experiments  at  low  velocities  in  Reference  4. 

Calculation  of  the  vane  response  characteristics  were  made  for  shock 
waves  of  103  kPa  (15  psi) ,  207  kPa  (30  psi),  and  345  kPa  (50  psi)  at  an 
ambient  pressure  of  86.2  kpa  (12.5  psi),  an  ambient  density  of  1.074  kg/m^, 
and  a  sound  velocity  of  335  m/s  (1100  ft/s).  The  calculations  were  made  for 
conditions  immediately  behind  the  shock  front,  so  that  density  and  flow 
velocities  were  obtained  from  shock  front  relations  given  in  Reference  6.  The 
results  are  given  in  Table  1. 

The  value  of  tg  listed  is  that  time  for  the  amplitude  of  the  oscillation 
to  decrease  by  one-half,  i.e..  In  Bq^^R  “  0.693.  The  response  times  are  long 
compared  to  the  usual  response  times  sought  for  blast  wave  measurements,  and 
are  longer  still  for  a  further  decrease  in  oscillation  amplitude. 

The  experimental  data  in  Reference  4  show  that  the  damping  for  low  speeds 
may  be  about  twice  as  large  as  that  predicted  by  equation  (12)  for  a  delta¬ 
shaped  vane.  However,  even  with  response  times  reduced  by  a  factor  of  two, 
they  are  longer  than  desired.  A  factor  that  will  complicate  interpretation  of 
a  record  from  such  a  vane  in  flow  changing  rapidly  in  flow  direction  is  that 
the  response  time  depends  on  the  velocity  of  the  flow. 

Another  measure  of  vane  performance  is  the  "direction  constant,"  defined 
in  Reference  7  as  the  length  of  air  passage  moving  at  speed  U  before  the  vane 
can  attain  one  cycle  in  its  direction  response.  It  is  given  by  the  relation 


T  U 
n 


21T  ( 


pAra'' 


(17) 


where 


T 

n 


1/f 

n 


2"  /  (») 

n 


(18) 
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Values  for  and  TnU  are  given  in  Table  1.  The  values  of  TnU  are 
approximately  one  metre,  but  they  are  not  constant.  Because  of  the  low 
damping  the  vane  will  require  several  cycles  to  settle,  so  the  length  of  air 
passage  at  speed  U  before  the  vane  is  near  equilibrium  will  be  about  three 
metres  or  more. 

Because  of  the  slow  response  time  and  low  spaclal  resolution  the  movable 
vane  sensor  of  the  type  shown  in  Figure  5  does  not  seem  to  be  suitable  for  use 
in  blast  waves. 
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Table  1.  Vane  response  characteristics  at  shock  front  conditions. 
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SECTION  5 


ALTERNATIVE  FLOW  DIRECTION  SENSORS 

A  technique  for  determining  flow  direction  and  magnitude  Is  that  of 
placing  a  body  of  simple  shape  In  the  flow  and  measuring  the  components  of 
force  acting  upon  It.  Reference  8  describes  a  small  rectangular  cantilever 
beam  anemometer.  It  Is  fragile  and  presents  a  changing  area  to  the  flow  as 
flow  angle  varies,  and  therefore  Is  unsuitable  for  the  present  application. 
Spheres  have  been  used  for  air  and  water  flow  measurements  (Ref.  9).  Spheres 
were  used  on  past  nuclear  tests  and  on  HE  tests  with  limited  success.  A 
cantilever  beam  gage  with  a  sensitive  center  section  was  used  on  HE  tests 
(Refs.  10  and  11).  The  response  times  of  these  gages  are  too  long  for  the 
application  of  Interest. 

A  cylindrical  cantilever  beam  mounted  on  a  single  component  balance  Is 
described  in  Reference  12.  Reference  13  describes  a  hollow  cylindrical 
cantilever  beam  designed  for  use  In  a  blast  wave.  Two  versions  were  built 
with  natural  frequencies  of  500  and  1000  Hz  and  tested  in  a  shock  tube,  where 
the  performance  was  judged  to  be  satisfactory. 

One  gage  was  0.0572  m  (2.35  in.)  in  diameter  and  0.328  m  (12.9  in.)  long, 
and  the  other  was  the  same  diameter  but  0.231  m  (9.10  in.)  long.  They  were 
Instrumented  for  longitudinal  strain  with  strain  gages  spaced  at  Intervals  of 
90  degrees  around  the  base  of  each  beam.  These  gages  were  large  and  slow,  but 
the  concept  seems  most  appropriate  of  those  found  in  the  literature  for  the 
current  application.  Therefore  a  small  hollow  cylindrical  cantilever  beam 
will  be  considered  as  a  means  of  Indicating  flow  direction. 


19 


SECTION  6 


CANTILEVER  BEAM  FLOW  DIRECTION  SENSOR 

The  proposed  configuration  for  using  the  hollow  cylindrical  beam  as  a 
flow  direction  sensor  Is  shown  In  Figures  6  and  7.  The  configuration  shown  on 
the  left  In  Figure  6  measures  flow  direction  in  the  horizontal  plane  and  In 
the  vertical  plane  passing  through  the  gage  and  ground  zero.  The  output  of 
the  beams  on  the  same  axis  are  added,  providing  higher  signal  levels  and  some 
compensation  for  the  lack  of  symmetry  and  partial  shielding  effects  for  off 
angle  flow.  Figure  7  shows  a  cross-section  of  this  configuration. 

To  measure  flow  direction  in  the  vertical  plane  perpendicular  to  the 
radius  from  ground  zero  will  require  a  beam  mounted  separately,  as  shown  to 
the  right  In  Figure  6.  The  end  loading  on  the  beam  will  be  significant,  and 
higher  sensitivity  to  bending  will  be  required.  A  beam  with  only  a  sensitive 
center  section  (Refs.  10  and  11)  may  be  required  for  this  orientation. 

Strain  gages  were  selected  as  the  means  of  measuring  beam  deflection 
because  they  represent  a  well  known  technology  and  can  be  connected  easily  so 
that  the  output  of  the  two  beams  on  the  same  axis  are  added  together.  The 
strain  gages  are  mounted  to  measure  longitudinal  strain  at  90  degree  Intervals 
around  the  inside  surface  of  the  beam  at  Its  base. 

The  beams  are  made  as  small  as  can  be  Instrumented  using  strain  gages. 
This  corresponds  to  a  diameter  of  about  0.010  m  (0.4  in.)  if  the  strain  gages 
are  mounted  on  the  outside  of  the  cylinder,  and  0.0127  m  (0.5  in.)  if  the 
strain  gages  are  mounted  inside  the  cylinder.  The  two  factors  affecting  size 
are  the  requirement  for  adequate  access  to  Install  the  gages  and  a  minimum 
curvature  of  the  surface  on  which  the  gages  are  to  be  mounted  (Ref.  14). 

The  output  from  each  beam  will  be  two  signals  indicating  the  magnitude  of 
bending  and  hence  force  In  two  perpendicular  directions,  so  that  force 
direction  and  magnitude  in  the  plane  perpendicular  to  the  axis  can  be 
determined.  Because  of  the  varying  drag  coefficient  across  the  subsonic. 
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transonic,  and  supersonic  flow  regimes  and  presence  of  dust  In  the  flow,  the 
force  component  measured  cannot  be  used  to  make  a  quantitative  determination 
of  dynamic  pressure. 


To  obtain  estimates  of  the  strain  at  the  base  of  the  beam,  simple  beam 
bending  relations  are  used  (Ref.  15). 

Strain  In  a  fiber  In  a  bent  beam  Is 

e  *  h/R  (19) 


where  h  =  distance  from  neutral  axis,  center  of  beam, 
R  =  radius  of  curvature  of  the  beam. 


Now 

1/R  =  (^)  /  (1  +  (£)^)^^^ 
dx^ 


(20) 


where  y  =  f(x),  the  function  describing  the  deflection  y  of  the  beam 
versus  x  distance  along  Its  length. 


The  slope  dy/dx  Is  very  small,  so  that  Its  square  In  equation  (20)  can  be 
neglected,  so  that 

1/R.  4 


dx 


(21) 


and 


E  =  h 


dx^ 


(22) 


For  a  uniformly  distributed  load  on  a  cantilever  beam  (Ref.  15), 


x^  (x^  -  4K.X  +  61^) 


(23) 
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where 


w  =  force  per  unit  length 

I  =  length  of  beam 

X  >  distance  from  supported  end. 


,2 

d  y  _  w 

,  2  ■  24EI 
dx 


(12x2  -  2hlx  +  121^) 


Then 

e  =  h  •  2!^  (12x^  -  24ix  +  12£^) 


For  X  =  0,  the  strain  at  the  base  is 


e 


wh 

2EI 


i 
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In  a  blast  wave  flow,  w  =  CpA’Q/i 

where  Cq  =  drag  coefficient 
A  =  area  =  J^^dg 
dg  =  outside  diameter  of  beam 
Q  =  dynamic  pressure  =  p\5^/2 

For  strain  gages  on  the  Inside  of  the  beam,  h  =  di/2,  where 
interior  diameter  of  the  beam. 

Substituting  these  values  in  equation  (26)  produces  the 
relation  for  the  maximum  strain  on  the  inside  wall: 


e 


=  s  ^ 


d  d 
o  1 

4  El 


(24) 


(25) 


(26) 


di  is  the 

following 


(27) 
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The  natural  frequency  of  vibration  of  a  cantilever  beam  of  uniform 
section  is  given  by  (Ref.  16) 


f  =  (3.52  /  2tt)  (El/mi^)^ 
n 


(28) 


where  m  =  mass  per  unit  length. 

For  the  hollow  cylindrical  beam  of  interest  here 


I  =  (tt/64)  (do^  -  di^) 


Let  Dq  =  dg/i  and  =  d^/i. 

Then  equation  (28)  can  be  written  as  follows: 

=  3-52  ^ 

n  8tt£.  P 

m 


<■>0*  - 
(d/  -  d/) 


h. 


(29) 


(30) 


where  Pjn  =  density  of  beam  material. 


The  proposed  material  for  the  beams  is  17-4PH  stainless  steel,  condition 
H900  (Ref.  14).  Using  E  =  200  x  10^  kPa  (29  x  10^  psi)  and  p^  =  7-835  x  103 
kg/m3  (0.283  lb/in.3),  equation  (30)  becomes  as  follows: 


f 

n 


707.5 

I 


(31) 


where  I  is  in  meters,  and 
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f 

n 


27860 

I 


(32) 


where  1  is  in  inches. 

The  beam  will  be  highly  underdamped,  so  that  the  usable  upper  frequency 
will  be  about  0.25  f^. 

The  output  signal  from  the  gage  should  be  filtered  prior  to  recording. 
The  gage  system  should  be  set  up  so  that  the  oscillations  remain  linear,  and 
the  electrical  output  remains  linear  and  it  is  not  clipped  prior  to  filtering. 

The  proposed  cantilever  beam  has  the  following  dimensions: 

i  =  0.0406  m  (1.60  in.) 

do  =  0.0147  m  (0.58  in.) 

di  =  0.0127  m  (0.50  in.) 

The  maximum  strain  on  the  inside  wall  from  equation  (27)  is  .000257.  For 
a  value  of  Q  of  689  kPa  (100  psi) ,  Cq  =  1.0,  E  =  200  x  10^  kPa  (29  x  10^  psi) , 
and  the  dimensions  listed  above,  equation  (27)  produces  a  strain  of  .000257. 

For  semiconductor  strain  gages  with  a  gage  factor  of  150  in  a  two-arm 
bridge  to  measure  bending,  the  output  from  a  beam  will  be  19.3  mv/v.  At  5 
volts  excitation  the  output  will  be  96.5  millivolts,  which  is  adequate. 

The  natural  frequency  of  the  beam,  calculated  from  equation  (31),  is  8330 
Hz.  Because  of  the  low  damping  the  upper  usable  frequency  will  be  about  0.25 
fjj  =  2080  Hz,  with  a  period  of  0.00048  s  and  a  rise  time  at  about  0.00012  s. 
This  time  is  short  enough  to  provide  a  useful  sensor  according  to  the 
discussion  in  Section  2  of  this  report. 
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The  cost  estimates  for  prototype  gages  are  as  follows: 


2-Beam 

4-Beam 

Hardware : 

Beams 

$  280 

$  560 

Beam  Probe 

$  530 

$  810 

$  540 

$1100 

Strain  Gaging 

$1000 

$2000 

[$500/beam] 

Static  Calibration 

$  300 

$  600 

TOTAL 

$2110 

$3700 

These  costs  apply  for  gage  designs  that  have  been  evaluated  and  found  to 
be  satisfactory.  No  development  costs  are  included. 

The  2-Beam  version  will  measure  flow  direction  in  the  vertical  plane  in 
the  line  through  the  gage  and  ground  zero.  The  4-Beam  version  will  measure 
flow  direction  in  the  horizontal  plane  and  in  the  same  vertical  plane  as  the 
2-Beam  version. 
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SECTION  7 


RECOMMENDED  TEST  PROGRAM 


Evaluation  and  certification  of  the  flow  direction  sensor  will  require 
static  testing,  and  tests  In  shock  tubes  and  wind  tunnels.  In  static  tests 
Its  output  versus  direction  of  force  applied  will  be  measured  and  Its  angular 
resolution,  symmetry  of  response,  linearity  and  hysteresis  determined.  After 
the  static  tests  have  been  completed  a  series  of  dynamic  tests  In  both  shock 
tubes  and  wind  tunnels  will  be  conducted. 

7.1  STATIC  CALIBRATION. 


For  the  static  calibration  a  concentrated  transverse  load  P  will  be 
imposed  on  the  end  of  the  cantilever  beam,  as  shown  in  Figure  8.  The  angle 
0  is  the  angle  between  the  line  of  action  of  the  load  P  and  the  strain  gage 
location.  For  a  selected  angle  0,  the  outputs  from  the  strain  gage  bridges 
will  be  recorded  as  the  load  is  increased  in  steps  from  zero  to  maximum  and 
decreased  in  steps  to  zero.  This  procedure  will  be  followed  for  different 
values  of  0.  The  data  obtained  will  be  examined  to  evaluate  the  linearity  and 
symmetry  of  strain  response  versus  angle  of  loading. 

The  data  will  also  be  compared  with  a  relation  derived  using  equation 
(22).  For  the  cantilever  beam  with  a  concentrated  transverse  load  P  at  its 
free  end,  the  deflection  is  (Ref.  15) 


y 


Px 

6EI 


(3)1  -  x^). 


(33) 


and 
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d  y  ^  ^ 

^  2  El 
dx 


(1 


x) . 
e 


(34) 


From  equation  (22)  the  bending  strain  is 
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Figure 


8.  Schematic  of  transverse  load  P  and  strain  gage  locations  for  static 
calibration. 
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e  =  h 


^2  d.P 

2  2EI  ^ 

dx 


(35) 


At  the  base  of  the  beam  and  in  terms  of  the  component  of  P  acting  at  the 
strain  gage  location,  the  bending  strain  is 


e  = 


dj^P  cos  6  I 

m 


(36) 


The  data  obtained  in  the  static  calibration  will  be  compared  with  values 
computed  using  equation  (36) . 

The  data  will  be  examined  to  determine  the  angular  resolution,  symmetry 
of  response,  linearity  and  hysteresis  of  the  gage. 

7.2  SHOCK  TUBE  TESTS. 

In  calibrating  the  flow  direction  gage  versus  incident  angle  of  flow, 
only  flows  incident  in  the  forward  hemisphere  aligned  with  the  gage  and 
support  shaft  are  considered.  Flows  not  incident  in  this  hemisphere  will  be 
obstructed  by  the  support  shaft  for  the  cantilever  beam. 


Figure  9  shows  the  proposed  incident  angles  in  the  vertical  plane 
containing  the  vertical  beam  axis  and  support  shaft.  The  angular  interval  is 
22.5  degrees,  and  the  total  number  of  incident  angles  for  testing  is  nine. 
Figure  10  shows  the  four  planes  in  which  the  incident  flow  is  to  be  varied. 
Mapping  of  gage  response  in  the  forward  hemisphere  requires  36  tests. 


The  first  test  series  is  to  be  conducted  in  the  0.6  m  air  driven  shock 
tube  at  the  Ballistic  Research  Laboratory  (BRL),  Aberdeen  F  ovlng  Ground, 
Maryland.  A  step  shock  of  172  kPa  (25  psi)  or  the  maximum  allowed  will  be 
used  for  all  tests.  The  dynamic  response  characteristics  of  the  gage  will  be 
determined,  and  filtering  requirements  investigated.  The  data  will  be 
analyzed  to  determine  angular  resolution  and  symmetry  of  response. 
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Figure  9.  Flow  angles  for  calibration  in  the  vertical  plane  containing  the 
vertical  beam  axis  and  probe  support. 
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Figure  10.  Planes  at  0,  45,  90,  and  135  degrees  in  which  flow  incident  angles 
are  to  ne  varied  for  calibration. 


After  these  initial  tests  and  the  wind  tunnel  tests  are  completed 
satisfactorily,  the  gage  will  be  installed  in  the  0.43  m  (17  in.)  shock  tube 
at  TRW  for  a  more  limited  test  series.  This  shock  tube  can  produce  higher 
shock  pressure  than  the  BRL  tube.  Tests  can  be  conducted  with  and  without 
dust  in  the  flow. 

The  Initial  test  in  the  TRW  tube  will  be  made  at  the  same  shock  pressure 
used  in  the  shock  tube  tests  in  the  BRL  tube  to  check  the  correlation  of  gage 
response.  Then  tests  will  be  conducted  at  345  kPa  (50  psi) ,  483  kPa  (70  psi) , 
and  621  kPa  (90  psi).  Four  shots  are  to  be  fired  at  each  pressure  level.  Two 
will  be  with  the  gage  aligned  with  the  flow  with  and  without  dust.  Two  will 
be  with  the  gage  at  an  off  angle  of  45  degrees,  with  and  without  dust.  The 
data  will  be  examined  for  differences  between  response  with  and  without  dust, 
and  between  pressure  levels  and  previous  test  results. 

The  design  of  the  TRW  test  series  may  be  changed  based  on  results  of  the 
BRL  tube  tests  and  wind  tunnel  tests. 

7.3  WIND  TUNNEL  TESTS. 

The  tests  in  the  BRL  shock  tube  described  above  will  provide  gage 
response  in  a  subsonic  flow.  Only  supersonic  wind  tunnel  testing  is  proposed 
because  of  the  high  velocity  of  the  flows  in  the  precursed  blast  wave  where 
measurement  of  flow  direction  is  of  interest. 

Tests  are  to  be  conducted  at  Mach  numbers  of  1.5,  2.0,  and  2.5  over  the 
full  range  of  incident  angles  in  the  forward  hemisphere  shown  in  Figures  9  and 
10.  The  data  will  be  examined  for  angular  resolution  and  symmetry,  and 
compared  with  the  data  from  the  tests  in  the  BRL  shock  tube. 
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SECTION  8 


CONCLUSIONS 


The  precursor  flow  field  contains  regions  with  large  and  fast  changes  In 
flow  direction  relative  to  a  fixed  sensor.  Examination  of  a  numerical 
calculation  of  a  precursed  blast  wave  Indicated  a  need  for  a  sensor  with  a 
response  time  of  0.2  ms  or  less. 

The  response  time  for  a  rotating  wind  vane  sensor  Is  Inversely 
proportional  to  the  flow  velocity.  For  what  was  considered  to  be  the  smallest 
practical  size  for  a  vane  (22  mm),  the  response  times  were  too  slow  for  use  In 
precursed  blast  waves. 

After  examining  alternative  means  of  sensing  flow  direction,  the 
direction  of  bending  of  a  cylindrical  hollow  cantilever  beam  protruding  Into 
the  flow  was  selected  as  the  best  means  for  Indicating  flow  direction.  It  was 
possible  to  design  a  beam  with  a  response  time  short  enough  for  use  In 
precursor  experiments.  The  beam  Is  a  hollow  cylinder  41  mm  long  and  15  mm  In 
diameter  with  a  natural  frequency  of  8300  Hz.  Because  It  Is  underdamped,  the 
upper  usable  frequency  Is  about  2080  Hz,  with  an  associated  rise  time  of  0.12 
ms . 


The  basic  cost  of  the  cantilever  beam  flow  direction  sensor,  once 
developed.  Is  estimated  to  be  $2100  for  sensing  In  one  plane,  and  $3700  for 
sensing  In  two  perpendicular  planes. 

For  development  and  certification  of  the  sensor,  static,  shock  tube,  and 
wind  tunnel  tests  will  be  required. 
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APPENDIX 
LIST  OF  SYMBOLS 


rate  of  change  of  lift  coefficient  with  respect  to  3  evaluated  at 

e  -  0 

projected  vane  area  in  the  plane  of  zero  angle  of  attack,  or  the 
projected  area  of  a  cantilever  bean 

aspect  ratio  of  vane 

drag  coefficient  for  cylindrical  cantilever  beam 
lift  coefficient  for  vane 

interior  diameter  of  hollow  cylindrical  beam 
outside  diameter  of  cylindrical  befun 
aerodynamic  damping  constant  for  vane 

ratio  of  interior  diameter  to  length  for  a  hollow  cylindrical  beam 
ratio  of  outside  diameter  to  length  for  a  hollow  cylindrical  beam 
Young's  Modulus  for  beam  material 
natural  frequency  of  oscillations 

aerodynamic  force  acting  normal  to  vane  at  its  center  of  pressure 
distance  from  neutral  axis  of  a  beam 
length  of  beam 

mass  per  unit  length  of  beam 
unit  angle  torgue  acting  on  vane 

concentrated  transverse  load  at  end  of  cantilever  beam 
dynamic  pressure  of  flowing  air 

distance  between  vane  center  of  pressure  and  the  rotation  axis 

radius  of  curvature  of  beam 

time 

vane  response  time,  the  time  for  the  amplitude  of  the  envelope  of 
oscillations  to  reduce  by  a  specified  amount 
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t{^  -  time  for  vane  oscillations  to  decrease  in  amplitude  by  one -half 

-  period  of  oscillations  at  the  natural  frequency 
U  -  wind  velocity  vector 

U  -  wind  velocity  magnitude 

Uy  -  air  velocity  relative  to  vane 

w  -  force  per  unit  length  along  a  beam 

X  distance  along  beam 

y  -  deflection  of  cantilever  beam 

g  -  angle  of  wind  vector  with  respect  to  vane 

gQ  ~  peak  amplitude  of  angular  deviation  of  vane 

g{(  -  vane  oscillation  amplitude  at  tp 

By  -  angle  between  air  velocity  vector  and  vane,  the  instantaneous 
angle  of  attack 

e  =  strain  in  a  fiber  of  a  bent  beam 

t,  ”  damping  ratio 

0  =  angle  describing  orientation  around  axis  of  cylindrical  cantilever 

beam 

p  =  air  density 

p,„  =  density  of  beam  material 

Wp  =  angular  natural  frequency 
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ATTN:  APALLONE 
ATTN:  WREINECKE 

BOM  CORPORATION 

ATTN;  LIBRARY 

BDM  INTERNATIONAL  INC 
ATTN:  EDORCHAK 

BOEING  CO 

ATTN:  ALUNDE 
ATTN:  RJWOOD 

CALIFORNIA  RESEARCH  4  TECHNOLOGY.  INC 
ATTN:  J  THOMPSON 

CARPENTER  RESEARCH  CORP 
ATTN:  HJ  CARPENTER 


Dist'2 


ODM-TR-M-IM  (DL  CONTNIUfD) 


DAYTON.  UNIVERSITY  OF 
ATTN:  BWILT 
ATTN:  N  OLSON 
ATTN:  R  SERVAIS 

GENERAL  RESEARCH  CORP 
ATTN:  RPARISSE 

GEO  CENTERS,  INC 

ATTN:  E  MARRAM 

GEORGIA  INSTITUTE  OF  TECHNOLOGY 
AHN:  C  BROWN 

H  &  H  CONSULTANTS.  INC 
ATTN:  WHALL 

IIT  RESEARCH  INSTITUTE 

ATTN:  DOCUMENTS  LIBRARY 

KAMAN  SCIENCES  CORP 
ATTN:  LMENTE 
ATTN:  LIBRARY 

KAMAN  SCIENCE^  CORP 

ATTN:  BKINSLOW 

KAMAN  SCIENCES  CORP 
ATTN:  E  CONRAD 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

KAMAN  SCIENCES  CORPORATION 
ATTN.  DASIAC 

KARAGOZIAN  AND  CASE 

ATTN:  J  KARAGOZIAN 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  L  CHASE 
ATTN:  R  SMITH 
AHN:  S  SALISBURY 
ATTN:  T  FISHER 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  TECH  INFO  CTR 

MISSION  RESEARCH  CORP 

ATTN:  TECH  LIBRARY 

NICHOLS  RESEARCH  CORP,  INC 
ATTN:  RBYRN 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  HBRODE 

PACIFICA  TECHNOLOGY 

ATTN:  TECH  LIBRARY 


R  &  D  ASSOCIATES 

ATTN:  CKBLEE 
ATTN:  D  SIMONS 
ATTN:  P  RAUSCH 
ATTN:  T  MAZZOLA 

ATTN:  TECHNICAL  INFORMATION  CENTER 

RAND  CORP 

ATTN:  BBENNEH 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  L  PINSON 

S-CUBED 

ATTN:  RDUFF 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  K  SITES 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  TECHNICAL  UBRARY 
ATTN:  W  PLOWS 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  MCRARY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  TECHNICAL  UBRARY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  COCKAYNE 
ATTN:  RSIEVERS 
ATTN:  WCHADSEY 
AHN:  WLAYSON 

SCIENCE  APPLICATIONS  INTL  CORP 
AHN:  GBINNINGER 

SRI  INTERNATIONAL 

ATTN:  DKEOUGH 
ATTN:  D  MCDANIEL 

TECH  REPS.  INC 
2CYS  ATTN:  G  WATSON 
2  CYS  ATTN:  N  ETHRIDGE 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  JRAVENSCRAFT 
ATTN:  P  SHELTON 

TRW  INC 

ATTN:  JTAMBE 
ATTN:  TECH  INFO  CTR 

TRW  SPACE  &  DEFENSE  SYSTEMS 
ATTN:  W  WAMPLER 


PHYSICS  APPLICATIONS.  INC 

ATTN:  DOCUMENT  CONTROL 


PHYSICS  INTERNATIONAL  CO 
ATTN:  JSHEA 


Dist-3 


